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The study of motions of small particles suspended in a fluid has always been an interesting 
subject in physics. 1 The dynamics of the particles is determined by the statistical properties of the 
random forces resulting from interactions between the particle and the surrounding fluid molecules. 
Brownian diffusion of small particles in a fluid at thermal equilibrium is one of the classical and 
best understood examples. Sedimentation, wherein heavy particles fall under the action of gravity 
through a fluid in which they are suspended, on the other hand, represents a “self-induced hy- 
drodynamic diffusion”. 2 The character of the sedimentation depends upon the interplay between 
long-range hydrodynamic forces, random Brownian forces and direct inter-particle forces. 3 The 
main objectives of our NASA research projects carried out at Oklahoma State University are: (1) 
study motions of colloidal particles under different random forces in the carrier fluid, and (2) un- 
derstand the effect of gravity on the particle motion when the density of the particles is different 
from the carrier fluid. 

In this paper we focus on the polymer-induced depletion attraction and its effect on colloidal 
sedimentation in colloid-polymer mixtures. We first report a small-angle neutron scattering (SANS) 
study of the depletion effect in a mixture of hard-sphere- like colloid and non- adsorbing polymer. 
Then we present results of our recent sedimentation measurements in the same colloid- polymer 
mixture. A key parameter in controlling the sedimentation of heavy colloidal particles is the inter- 
particle potential U(r), which is the work required to bring two colloidal particles from infinity 
to a distance r under a given solvent condition. This potential is known to affect the average 
settling velocity of the particles, 3 and experimentally one needs to have a way to continuously vary 
U(r ) in order to test the theory. The interaction potential U(r) cam be altered by adding polymer 
molecules into the colloidal suspension. In a mixture of colloid and non-adsorbing polymer, the 
potential U(r) can develop an attractive well because of the depletion effect 4 , in that the polymer 
chains are expelled from the region between two colloidal particles when their surface separation 
becomes smaller than the size of the polymer chains. The exclusion of polymer molecules from the 
space between the colloidal particles leads to an unbalanced osmotic pressure difference pushing 
the colloidal particles together, which results in an effective attraction between the two colloidal 
particles. The polymer-induced depletion attraction controls the phase stability of many colloid- 
polymer mixtures, which are directly of interest to industries. 

It has been shown that the potential U(r) has the form 4 
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where cf is the particle diameter, lip is the osmotic pressure of the polymer molecules, and R g is 
their radius of gyration. The volume of the overlapping depletion zones between the two colloidal 
particles is given by 4 



where v p — (4n/3 )R 3 g is the volume occupied by a polymer chain and A = 1 + 2R g /a. In th 
experimet, we measure the colloidal (partial) structure factor, S C (Q), which is directly related to th 
interaction potential U(r). An advantage of using SANS is that one can eliminate the undesirabl 
scattering from the polymer chains by using isotopically mixed solvents. (In the discussion below 
the subscripts c and p will be used to refer to colloid and polymer, respectively.) 

The colloidal particles used in the experiment consisted of a calcium carbonate (CaCOz) cor 
with an adsorbed monolayer of a randomly branched calcium alkylbenzene sulphonate surfactan 
These particles have been well characterized previously using SANS and small- angle X-ray seal 
tering (SAXS) techniques. 5 Our recent SANS and SAXS measurements revealed that the particl 
has a core radius, Ro = 2.0 nm, and a monolayer thickness, 6 = 2.0 nm. Previous dynamic ligl 
scattering experiment 6 has shown that the colloidal particles are relatively monodispersed wit 
~ 10% standard deviation in particle radius. The polymer used in the study was hydrogenate 
polyisoprene (poly-ethylene-propylene or PEP), a stable model polymer. The molecular weigl 
of the PEP was M p = 26,000. Decane has been found to be a good solvent for both the colloi 
and PEP. 6 Because decane and PEP are both protonated, the polymer chains in the mixture ai 
invisible to neutrons. Our Zimm analysis of the SANS data from the pure PEP/deuterated-decar 
solution has shown that the polymer chains have a radius of gyration R g = 8.3 nm and their secon 
virial coefficient A 2 M P = 44.4 (cm* / gm). With the measured A 2 one can define an effective har( 
sphere radius Rh s via 4(4n/3)R* hs = A 2 M Thus we have Rh s = 4.8 nm, which agrees well wit 
our previous light scattering measurement. 6 The light scattering experiment has revealed that tl: 
PEP chains do not adsorb onto the colloidal surfaces, and the phase separation in the colloid-PE 
mixture samples occurs at the concentrations very close to the depletion prediction. 6 Because tl 
basic molecular interactions are tuned to be simple, the SANS measurements in the colloid-PE 
mixture can be used to' critically examine the current depletion theory. The SANS measuremen 
were performed at the High Flux Beam Reactor in the Brookhaven National Laboratory. The ii 
cident neutron wavelength A 0 = 7.05 ± 0.4 A, and the usable range of the scattering wave numb< 
Q [= (47r/Ao)sin(0/2), with 8 being the scattering angle] was 0.007 A" 1 < Q < 0.15 A -1 . Tl 
structure factor S C (Q ) was obtained using the equation S C (Q) = I(Q)/[pcPc(Q)\ , where I(Q ) 
the scattered intensity of the mixture samples, p c is the colloid number density, and P C (Q) is tl 
scattering intensity per unit concentration measured in a dilute pure colloidal suspension, in whic 
S C (Q) = 1. All the scattering measurements were conducted at room temperature. 

To reduce the fitting ambiguity and pinpoint the control parameters for the depletion effec 
we prepared three series of mixture samples with the colloid volume fraction <f) c — 0.146, 0.0£ 
and 0.038, respectively. For each series of the samples, <f> c was kept the same and the polyun 
concentration C p ( gm/cm *) was increased until the mixture became phase separated (except fc 
the series with <j> c = 0.038) with a visible interface, which separates the dark brown colloid-ri< 
phase from the light brown colloid-poor phase. Fig. 1 compares the measured S C (Q) for thn 
values of C p when (a) <j) c = 0.146 and (b) (f> c = 0.086. It is seen that the main effect of adding PE 
into the colloidal suspension is to increase the value of S C (Q) in the small-Q region, whereas tl 
large-Q behavior of S C (Q) remains nearly unchanged. The solid curves in Fig. 1 show the calculate 
S C (Q) using U(r) in Eq. (1). Note that there are four fitting parameters in the calculation: tl 
diameter < 7 , the volume fraction 4> c , the dimensionless interaction amplitude P — U p v p / J cbT, ai 
the range parameter A. The first two parameters are used to describe the hard core part of U(r 
and the last two parameters are used for the attractive tail. It is found from the fitting that f< 
a fixed colloid concentration, the fitted values of a and <j> r for the mixture samples do not chanj 
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very much with C p , and they are very close to those obtained from the corresponding pure colloidal 
suspensions. Furthermore, the fitted A also remains constant for different <j> c and C p , and its best 
fit value is A = 2.9. This value is close to the calculated A = 1 + R g /(Ro + 6) = 3.07. With 
the above three fitting parameters fixed, we were able to fit all the scattering data from different 
mixture samples (19 samples in total) with only one free parameter - the interaction amplitude P. 

Figure 2a shows the fitted P as a function of the effective polymer volume fraction <f> p = 
C p /C * , where C* = M p /[(4ir /tyR^] is the polymer overlap concentration. It is seen that P first 

increases linearly with <j) p up to <f> p ~ 1 and then it levels off. For a given <f> p , P also depends 
upon 4>c- If the polymer molecules in the mixture are treated as an ideal gas, their osmotic 
pressure lip = ripkeT and hence P — U p v p /(kBT) = <f> p . Recently, Lekkerkerker et al. 7 pointed 
out that the polymer number density n p should be defined as n p = N p /V /, where N p is the 
total number of the polymer molecules and Vf = ot{4> c )V- is the free volume not occupied by the 
colloidal particles and their surrounding depletion zones. They have calculated a(<f> c ) as a function 
of <j) c . It is seen from Fig. 2b that once <f> p is scaled by the calculated a(<f > c ), 8 the three curves 
in Fig. 2a collapse into a single master curve. The solid curve in Fig. 2b is the fitted function 
P = -0.054 + 0.178(<^p/a) - 0.0245(<£ p /a) 2 . 

The fitted P consists of three terms. The small negative intercept indicates that there is a 
weak repulsive interaction between the soft surfactant shells of the colloidal particles. (To have a 
meaningful comparison with the fitted P for the mixture samples, we used the same U(r) as in Eq. 
(1) but change the sign of U(r) for r > a to fit the repulsive tail for the pure colloidal samples.) 
The linear coefficient should be unity for non-interacting polymer chains (an ideal gas), but our 
fitted value is 0.178. One plausible reason for the deviation is that with the effective potential 
approach, the polymer molecules axe assumed to be smaller than the colloidal particles and their 
number density should be much higher than that of the colloidal particles. In our experiment, 
however, these two assumptions are not strictly satisfied, and thereby the overlap volume Vo(r) in 
Eq. (2) is over-estimated. As a result, the fitted P becomes smaller than its actual magnitude, 
because U(r) in Eq. (1) is proportional to the product of Vo( r ) and P. Another possibility is 
that in calculating <f> p , a smaller characteristic length than R g should be used for the polymer 
chains. For example, if Rh, is used to compute <j> p , the linear coefficient will be increased by a 
factor of ( Rg/Rhs ) 3 — 5.2. The polymer-polymer interaction, which gives rise to the quadratic 
term in the fitted P, can have two competing effects on the depletion attraction. It may either 
increase the osmotic pressure (and hence P) because the polymer chains have a positive second 
virial coefficient, or reduce the depletion attraction because it requires the system to do more 
work to expel the polymer molecules from the depletion zones. Fig. 2b clearly shows that the 
polymer-polymer interaction tends to reduce the depletion attraction. Similar suppression effects 
are also found in recent theoretical calculations of the depletion attraction between two parallel 
plates immersed in an interacting polymer (or particle) solution. 9 

The above SANS measurements clearly demonstrate the effectiveness of using a non-adsorbing 
polymer to control the magnitude as well as the range of the interaction between the colloidal 
particles. It is shown that the amplitude of the potential U(r) is controlled by the polymer con- 
centration C p , and the range of U(r) is proportional to the radius of gyration R g of the polymer 
chains. Because the depletion potential U(r) can be continuously varied in the experiment, the 
colloid-PEP mixture becomes an ideal system for the study of the interaction effect on the colloidal 
sedimentation. In the experiment to be described below, we measure the average settling velocity 
v c (C p ) of the colloidal particles as a function of C p . The addition of the polymer molecules into 
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the colloidal suspension can have two competing effects on v c (C p ). It can either reduce v c (C p ) 
because the viscosity of the mixture solution is increased, or increase v c (C p ) because the depletion 
attraction between the particles is increased. Experimentally, one can separate the two effects by 
changing the colloid concentration. For sufficiently dilute colloidal mixtures, the distance between 
the colloidal particles is so laxge that their mutual interaction can be ignored. In this case, adding 
polymer into the colloidal suspension only affects the viscosity of the solution. The effect of U(r) 
on v c (C p ) can be studied in the concentrated colloidal mixtures, once the sedimentation of the 
individual particles through a polymer solution is understood. 

We now discuss the measurements of v c (C p ) in the dilute colloidal mixtures. The CaCO\ 
particles have a hydrodynamic radius, Rh = 5.0 nm, and their density d c ~ 2.0 gm/cm 3 . The 
density of the PEP is d p = 0.856 gm/cm 3 . The solvent (decane) density is d 9 = 0.73 gm/cm 3 ant 
its viscosity rj o = 0.84 cp (at T = 25° C). The Stokes velocity of the particle, determined by the 
balance between the accelerational force and the viscous drag, is 


Vc(C p ) 


2 R 2 h (d c - d 3 )A 

c(Cp) 


(3 


where A is the acceleration and rj c (C p ) is the viscosity experienced by the colloidal particles in th< 
polymer solution at the concentration C p . Because the particles are very small, their sedimentatior 
under earth gravity (A = g) is unobservable. To increase the settling velocity of the particles, w( 
used a commercial ultracentrifuge, Beckman Model L8-70M. The sample cells were made of Ultra 
Clear solid polymer and their size (diameter x height) was 14 x 95 mm. The distance betweei 
the middle of the sample cell and the center of rotation was r = 11.3 cm. All the samples wen 
centrifuged at the rotation speed / = 35, 000 rpm for 4 to 6 hours depending on the sample viscosity 
The corresponding centripetal acceleration A = (27r f) 2 f was 1.5 x 10 8 cm/s 2 (1.5 x 10 5 #), which wa 
large enough to cause the colloidal particles to settle 1 to 6 cm towards the bottom of the cell. Afte: 
the centrifugation, a clear interface could be observed by eye in the initially uniform solution. Thi 
interface separated the upper clear solvent region from the lower dark-brown colloid-rich region 
The travelling distance h of the interface was measured by a low- magnification microscope mountec 
on a translational stage controlled by a micrometer. The settling velocity was computed from tin 
measured h via v c (C p ) = h/t, where t is the running time. All the measurements were conducte< 
at 22°C . Because the polymer density is very close to that of the solvent, the sedimentation o 
the polymer molecules was negligible. Therefore, the colloidal particles settled through a uniforn 
PEP/decane solution at rest. To reduce systematic errors in the experiment, we present th* 
sedimentation data in terms of the velocity ratio R c = v c (C p — 0 )/v c (C p ). As shown in Eq. (3) 
R c — rj c (Cp)/r}o for the dilute colloidal mixtures. 

Figure 3 shows the measured R c as a function of C p at <j> c = 0.014. The molecular weight o 
the PEP was M p = 17, 500 and its R g ~ 5.8 nm. 10 It is seen that the measured R c (C p ) first increase 
linearly with C p up to C p ~ 0.07 gm/cm 3 and then it turns up sharply. Similar behavior was als« 
observed for the measured R c (C p ) with other molecular weights. The crossover concentration C 
is found to be independent of M p . With this value of we estimate the correlation length (o 
the mash size) of the polymer solution to be £ ~ R g (C p /C *)~~ 3y/4 ~ 3.5 nm, which is close to th 
size of the colloidal particles. The solid curve in Fig. 3 shows the macroscopic viscosity of th 
same polymer solution measured independently by Davison et al. with a capillary viscometer. 1 
The functional form of the solid curve is i ] p /t]q = 1 + [^]C P + kn([r]]C p ) 2 , with the coefficient 
[rj] = 25 cm 3 /gm and ku = 0.335. The dashed line is the linear plot rj p /rio = 1 + [r)]C p . Fig. 3 thu 
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reveals that the colloidal particles in the polymer solution feel the single-chain viscosity when their 
size is smaller than f. The particles experience the macroscopic viscosity of the polymer solution 
when their size become much larger than £. 

We now discuss the effect of the colloidal interaction on R c (C p ). Fig. 4 shows the measured 
R c vs. C p for different <f> c . It is seen that for the two lowest colloid concentrations ( <f> = 0.014 
and 0.02), the measured f? c (C p )’s superimpose with each other, indicating that the interaction 
between the colloidal particles is negligible and thus the velocity ratio R c equals the viscosity 
ratio r] c (C p )/T]o . As <p c increases, the colloidal particles feel more and more depletion attraction 
and, therefore, their settling velocity at a fixed C p is increased ( R c is reduced). For a fixed 4> c , 
the depletion attraction between the particles also increases with increasing C p . Consequently, 
the colloidal settling velocity is increased. This is seen in Fig. 4 once the effect of increasing 
viscosity is divided out from the measured R c vs. C p . For the first time, the experiment provides a 
quantitative relationship between the settling velocity v c (C p ) and the inter-particle potential U(r). 
Further theoretical calculation of v c (C p ) using the depletion potential U(r) is underway. The above 
measurements reveal a great potential of using sedimentation to measure the interaction between 
the colloidal particles and other macromolecules in complex fluids. 
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Fig. 1 - Measured S C (Q) of the colloid-PEP 
mixtures for different C p when (a) <j> c = 0.146 and (b) 
= 0.086. The values of C p ( gm/cm 3 ) in (a) are: 
0.0039 (circles), 0.0165 (triangles), 0.0308 (squares), 
and those in (b) are: 0.0038 (circles), 0.0233 (trian- 
gles) and 0.0652 (squares). The solid curves are the 
calculated S C (Q) using U{r) in Eq. (1). 
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Fig. 2 - Fitted P as a function of (a) <f> p = 
Cp/C * and (b) 4> p /ct. The values of (f> c in the mixture 
samples are: 0.146 (circles), 0.086 (triangles) and 0.038 
(squares). The solid curve in (b) is the fitted function 
P = -0.054 + 0.178(<£p/a) - 0.0245(</> p /a) 2 . 



Fig. 3 -Measured R c as a function of C p at <j> c = Fig. 4 - Measured R c vs. C p for different <j> c . 

0.014. The solid curve is r) p /r} 0 = 1 + fa] C p + &//([??] C p ) 2 The values of <j> c are: 0.014 (triangles), 0.02 (circles), 

with [tj] = 25 cm 3 /gm and it// = 0.335. The dashed 0.086 (diamonds), and 0.146 (squares), 

line is the linear plot r} p /r? 0 = 1 + [*]] C p . 
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